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Abstract 

We study the stability of the Harrison-Perkins-Scott (HPS) mixing pattern, assumed to hold at 
some high energy scale, against super symmetric radiative corrections. We work in the framework 
of a reference minimal supergravity model (mSUGRA) where supersymmetry breaking is universal 
and flavor-blind at unification. The radiative corrections considered include both RGE running 
as well as threshold effects. We find that in this case the solar mixing angle can only increase 
with respect to the HPS reference value, while the atmospheric and reactor mixing angles remain 
essentially stable. Deviations from the solar angle HPS prediction towards lower values would 
signal novel contributions from physics beyond the simplest mSUGRA model. 
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I. INTRODUCTION 



The discovery of neutrino oscillations has indicated a very peculiar structure 

of lepton mixing [ol, quite distinct from that of quarks. These data have triggered a rush of 
papers attempting to understand the values of the leptonic mixing angles from underlying 
symmetries at a fundamental level. An attractive possibility is that the observed pattern 
of lepton mixing results from some kind of flavour symmetry, such as A4, valid at a some 
superhigh energy scale where the dimension-five neutrino mass operator arises 0]. 

Here we reconsider the Harrison-Perkins-Scott (HPS) mixing pattern ^ within a simple 
reference model approach. Our only assumption is that at the high energy scale the tree-level 
neutrino mass matrix m^J'^'^ is diagonalized by the so-called HPS matrix, taken as. 



HPS 



/ ^/2/3 1/VS \ 
\-l/V6 1/V3 I/V2 J 



(1) 



which corresponds to the following mixing angle values: 



tan ^ATM 


= tan^ 


= 1, 


sin^ 6'chooz 


= sm^e% 


= 0, 


tan^ 6'soL 


= tan^ e^^2 


= 0.5 



(2) 



These predictions which hold at high energies may be regarded as a good first approximation 
to the observed values 0] indicated by oscillation experiments y, 3, 0, 0|. The diagonal 
neutrino mass matrix can be written as rh^J'^'^ = f/^ps ' ^^J"^"^ ' ^hps = diag(mi, 1712, ms), so 



that the tree-level neutrino mass matrix becomes 

^ |mi + |m2 "l"^! + |"^2 + |"^2 ^ 

\-imi + im2 |mi + |m2 - |mi + |m2 + / 



m 



tree 



(3) 



This form corresponds to a specific structure for the dimension-five lepton number violating 
operator. For example, it constitutes the most general ansatz that follows from a basic 
A4 symmetry for the neutrino mass matrix and the quark mixing matrix [7]. One of the 
central open questions in neutrino physics is to identify the exact mechanism of producing 
Fig. [H As a first step, here we will adopt a model-independent approach of considering the 
implications of Eq. assuming only the evolution expected in flavor-blind softly broken 
minimal supergravity at unification. This will provide us with a reference value that can be 
useful in the future for treating different models of neutrino mass ^. 
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Figure 1: Dimension five operator responsible for neutrino mass. 



II. RADIATIVE CORRECTIONS 



It has already been noted that radiative corrections present in the Standard Model renor- 
malization group equations (RGEs), leave the HPS "reference" predictions essentially sta- 
ble In addition to Minimal Supersymmetric Standard Model RGE evolution, here we 
consider also the effect of one-loop threshold effects [ll]. We will first consider the evolution 
of the neutrino oscillation parameters that follow from Eq. ([3]), which covers both the cases 
of degenerate as well as hierarchical neutrino masses. The radiatively corrected neutrino 
mass matrix in this case becomes 



l-loop 



m 



tree 



+ 6' 



m 



tree 



(4) 



where 



/5, 

S. 



ee ^/^e 



\5eT Kt/ 



(5) 



The diagonal elements include the threshold correction and the RGE running 



where the RGE running effect is 12| 



"'a 
16772 



In 



M, 



GUT 



Ml 



EWSB 



(6) 



(7) 



In order to get the analytic expressions for the threshold corrections, we proceed as in 
Ref. fu^. However, now we do not neglect Yukawa couplings, taking into account the fact 
that right- and left-handed charged sleptons mix. Therefore, the analytic expressions for 
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the deltas are 

6 2 

^ (^TT* h TT* \(^TT. D^* U^TT. l, 

HI3+3) 



^al^ ~ lQ:^^9UAlRia - haUX2Ria+3){9UAlRip " h/3UA2Ri* 

i=l A=l 

X Bi{ml+,ml) , 

i=l A=l ^ ' 

2 2 



Xb' 



E E E - /^«f/i24.+3)^?f^Aiiv^B2ri?5aoK_,, . 

i=l A=l B=l 

= EEE^I^Ar^2-^?'Ar^ir|Ar^4pi?Li?gC^oo(^^^ , 

i=l A=l B=l 

where we have evaluated the Feynman diagrams at zero external momentum, which is an 
excellent approximation as the neutrino masses are tiny. Here S^^^^f^^ ,{oi,P = e,/i, r), are 
the contributions from the chargino/charged slepton diagrams in Fig. [2] (a,c), respectively, 
while are the contributions from the neutralino/sneutrino diagrams. The values of 

the Saf3^s, in Eqs. ([5]) and are the sum of the four contributions given above. Analogous 
contributions exist corresponding to the symmetrized terms in Eq. (jl]) , required by the Pauli 
principle, as displayed in Fig. [2] (b,d). In the above formulas, U and V are the chargino 
mixing matrices and ^^+, = 1, 2), are chargino masses, while N is the neutralino mixing 
matrix with rn^()^, (A = 1, ..,4), denoting the neutralino masses. Finally, the matrices R^^''' 
denote the slepton/sneutrino mixing matrices, respectively. The coupling constant of the 
SU (2) gauge group is denoted g and that of U{1) is g'. Here ha is the charged lepton Yukawa 
coupling in the basis where the charged lepton masses are diagonal. Furthermore Bi and 
Coo are Passarino-Veltman functions given by 



° 2 ' 2 VmIwsb/ Mt-iy ^ 



where t = m\/rn^ and 



2 2 2N_l/o,oAX "^0 ^ , -2r2(r2-l)ln(ri) + 2r2(ri-l)ln(r2) 



CooK, mf, mi) = -(3 + 2A,) - - In + 



4 VMIwsb; 8(ri-l)(r2-l)(ri-r2) 

(10) 

where ri = m\/m^ and r2 = nii/rn^. We have used dimensional regularization, with e = 4— n 
and n is the number of space-time dimensions. The term = (2/e) — 7 + 41n(47r), where 
7 is Euler's constant, is divergent as e ^ 0. 
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(a) 



(b) 



(c) 



X+/X° X+/X° 



X+/X° X+/X° 



Figure 2: Feynman diagrams responsible for neutrino mass radiative corrections. The blob indicates 
an effective Lagrangian term obtained from integrating out the heavy right-handed neutrinos 

III. CORRECTIONS TO MIXING ANGLES: NUMERICAL RESULTS 

We now describe our numerical procedure. In order to compute the magnitude of the 
radiative corrections expected in the HPS anzatz we work in the framework of a refer- 
ence minimal supergravity model approach, with universal flavor-blind soft supersymmetry 
breaking terms at unification. Therefore the off-diagonal elements in the matrix in Eq. ([5]) 
are all zero ^ 

4m = 4r = = ^i,e = ^re = ^rfi = ^ ■ (H) 

We first have used the SPheno package [l^ to calculate spectra and mixing matrices 
within mSUGRA within the ranges: M1/2, Mq, Aq G [100, 1000] GeV, Aq with both signs, 
tan/5 G [2.5, 50] and yU with both signs. Then we have calculated the RGE running, Eq. ([7j), 
and the threshold corrections, Eqs. (jS]). We have explicitly checked that the dominant 
contribution to 5^^, defined in Eq. ([6]), always comes from the threshold corrections for 
a = e, /i. Also for a = r, threshold corrections are usually more important than RGE 
running contributions, typically 

~O(10(-^'-3)), V« (12) 
^ Nonzero off-diagonal elements may arise due to running, see discussion. 
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while 



5e ~ C(10(-^l'-9)) 



6, ~ O(10(-^'-^)) 



6r ~ O(10(-^'-2)) . 



(13) 



Note that only for very large values of tan/3, the RGE effect 6r is slightly larger than the 
threshold corrrections Srr- Using these radiative corrections we have computed the delta 
matrix in Eq. ([5]) and inserted it in the neutrino mass matrix at 1-loop given in Eq. (jl]). We 
have then numerically diagonalized the 1-loop neutrino mass matrix in Eq. (jl]) in order to 
obtain the neutrino masses and mixing angles. 

Notice that the HPS scheme only fixes neutrino mixing angles. Thus, the neutrino masses 
are free parameters. Making use of this freedom, we have used an iterative procedure in order 
to choose the parameters mi, m2 and m^, so that the numerically calculated 1-loop neutrino 
masses are such that the solar and atmospheric squared-mass splittings Am^QL ^^Itm 
reproduce the current best fit point value. In our numerical calculation we concentrate on 
normal hierarchy. We will comment on the case of inverse hierarchy at the end of the next 
section. 

The numerically calculated atmospheric and reactor neutrino mixing angles at low ener- 
gies do not deviate significantly from its HPS reference value at high energies. Indeed, the 
numerical results are: 



However, the solar neutrino mixing angle can be significantly affected. In Fig. [31 we 
have plotted the maximum deviation of the solar angle from the HPS reference value for 
tan/5 G [2.5, 50], as a function of rriu-^, for both extreme CP parity combinations for m,^^ and 
mu2- same sign (left panel) and opposite sign (right panel). All the other CP possiblities 
lie in between these two extreme cases. As can be seen, the solar mixing angle remains 
essentially stable in the case of opposite CP signs, while deviations are maximal in the case 
of same CP signs. In this case, the solar mixing angle always increases with respect to the 
HPS value, irrespective of mSUGRA parameters. Moreover we can get a rough upper bound 
on m^, of order 



for the mSUGRA parameter values: M1/2 = 100 GeV, Mq = -Ao = 10^ GeV, /i > and 
tan/? = 2.5. Note that the upper bound is sensitive to the values of tan/5. For higher values 
of tan/5 the radiative corrections are larger, implying a more stringent bound on m^-^, as 
indicated by the upper boundary of the red (dark) band of the left panel in Fig. [31 Here 
we have fixed solar and atmospheric mass squared splittings at their best-fit values from 



tan2^^^™<tan2^°3 + O(10^^), 
sin2^chooz<sin2^?3 + O(10-^). 



(14) 




(15) 
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Ref. [g]. However, we have explicitly checked that the effect of letting Aml^^ and Amg^^ 
vary within their current 3a allowed range is negligible, i. e. the bands obatined at the 
extreme values almost coincide with the ones in Fig. [31 




[eV] 



m [eV] 



Figure 3: Upper bound for the solar mixing parameter tan^^soL, as a function of m^-^ (in eV), 
for tan/3 = 2.5 (lower border of the red band) and tan/3 = 50 (upper border of the red band). 
On the left panel, m,^-^ and have the same CP sign. On the right panel, mj^^ and rrij^^ have 
opposite CP sign. The neutrino mass splittings are assumed to have their best fit value from Q]. 
The horizontal band corresponds to the 3ct allowed range for tan^ ^sql Q] • 



IV. ANALYTICAL UNDERSTANDING 

The numerical results presented above can be understood analytically as follows. If we 
perform the original HPS rotation to the 1-loop neutrino mass matrix in Eq. |H we get: 



where 



- l-loop _ rrT _ ™ 1-loop _ tt 
''V — '-^HPS ''V "^l 



HPS 



f (1 + 5ii)mi 5;"2'^i + ^T2^2 S'^^'mi + ^^^^mgX 

^S^mi + 6T2'm2 (1 + 622)7712 (52l'm2 + ^sT^s 
\5Ts'7n, + 5Ts'ms 5Ts'm2 + 5^,'ms {1 + 5 33) 7ns ) 



(16) 



(17) 



"11 

S22 
S33 



Tfl } 
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273' 
1 



The matrix in Eq. ( |T71) should be nearly diagonal and its off-diagonal elements determine 
the deviations from tri-bimaximality. We define the following parameters characterizing the 
deviations from tri-bimaximality: 

o tan(2ejj) ^ i-ioop\ / ^ i-ioop\ ' 

(19) 

so that 



(20) 



^ATM 


= ^23 


— "^23 


'chooz 


= 6*13 


— f^l3 


6'sOL 


= Q\2 


— ^^12 



Substituting the matrix elements in Eq. ( |T71) into Eq. ( |T9|) . we get: 



(-1 - (522)^2 + (1 + (533)?7l3 ' 

^ ^r3^mi+^r3^m3 



;-l - 5ii)mi + (1 + 533)"^3 

^r2^mi + (5^2^7712 
(-1 - 5ll)mi + (1 + 522)'"^2 
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Taking into account that for mSUGRA the off-diagonal elements in the matrix in Eq. ([5]) 
are all zero, see Eq. flTT]) . the 5's in Eq. f[T5]) become 

5- = = = ^(25L - - ^;.) , ^^^^ 

The deviations of the neutrino mixing angles from the HPS value given in Eqs. (I2m23p then 
become 

^ '^23(^^2 + ^3) , . 

(-l-5°2)m2 + (l + 53°3)m3 ' ^'^^ 

(-l-5?i)mi + (l + 5ym3 ' ^ ^ 

If ei2, given in Eq. (1271) . is always positive, 6'sol can only increase, see Eq. (120|) . The denom- 
inator in Eq. ( 1271) can be approximated to 

(-1 - 5?i)mi + (1 + 5^2))^2 ^ -mi + m2 > (28) 

and hence, by assumption, is always positive. The sign of ei2 will be the sign of 5^2 given by 
Eq. flMl) . Considering the expressions for the deltas given in Eq. ([8]) and bearing in mind 
that the Passarino-Veltmann functions depend rather smoothly on their arguments, we can 
take them out of the sum. The following very rough estimations of the threshold corrections 
result 

5aa ^ ^{?>g\Br + 4Coo) + g'\B^ + 4Coo)) , (« = e, /x) , 

Y (29) 



327r2' 



5rr ^ —^{3g\B, + 4Coo) + g"{B, + 4Coo) + 2hlB 



where we have neglected the charged lepton Yukawa couplings for a = e, /i. Using 

lim ^ , \ , =-2, (30) 
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Eq. (129|) becomes 



5a 

6^ 



327r2 
327r^ 



(a = e,/i) 



(31) 



Therefore, the contribution of the threshold corrections to is roughly 



25. 



5^^l - (^T 



167r 



2 r 



(32) 



Besides the threshold correction contributions, one has also to consider the RGE running 
contribution. Here the dominant part obviously is 5r-, given in Eq. ([7]). The approximated 
expression for (5°2, defined in Eq. ( 127|) . is then 

1 hi \ ^ , f Mr.^^T\ 



AO ^ 
"12 — 



3^2 



(25, 



3v^l67r2 



-5i + In 



^GUT 
^EWSB / 



(33) 



Considering that in the limit where the slepton mass goes to infinity, the Passarino-Veltman 
function Bi behaves as 

,2 



hm B,{ml ,mn^-ln 

o Xa' Li-' o 



(34) 



XA 



one obtains, from Eq. (133|) . 
1 



5° 



12 



3^2 



hi 



3V2l67r2 



In 



GUT 



Ml 



EWSB 



-In 



m-j 



m 



XA 



(35) 



which is always positive, thus explaining why ei2 > 0. Note that although the threshold 
corrections are in general larger than the RGE contributions, in there is a cancellation 
among the threshold corrections so that the 5^- RGE contribution becomes the relevant term. 
We have numerically checked that 



c(io(- 



-6-3)N 



(36) 



This cancellation among the threshold corrections is the reason why the solar neutrino 
mixing angle can only increase with respect its HPS reference value. 

We now turn to the other two neutrino mixing angles. In the mSUGRA framework 
the deviations from the HPS predictions are much smaller than found for the solar mixing 
parameter, and fit within their current experimental 3cr allowed range given in Ref. [gI for 
acceptable rriy^ values. The reason for this can be understood from Eqs. fl25H27l) . On the 
one hand, the deltas on the numerators, given by Eq. ( 12^ . are all of the same order. For 
small values of rriy^ the deviations are all negligible, since they are all proportional to the 
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previous deltas. For large m,y-^ values the neutrino masses are nearly degenerate so that the 
numerators in Eqs. fl25H27l) are all of the same order. The denominators in Eqs. fl25H27p can 
be approximated as 

(-l-<^22)"^2 + (l + 533)m3 ^ms-ma , (37) 
(-1 - 5°i)mi + (1 + (5°3)m3 ~ ms - mi , (38) 
(-1 - 5?i)mi + (1 + 5^2)^2 ~ ma - mi . (39) 

Although these mass differences are very small, m3 — ma and — mi are larger than 
ma — mi, thus making ea3 and €13 smaller than cia- 

We now comment briefly on inverse hierarchy. As can be seen from Eqs. fl37tl39p . for 
inverse hierarchy, ma — m-i is still much smaller than m3 — ma or m3 — mi, while the latter 
two just change sign but not the magnitude. We therefore expect that the above discussion 
remains essentially correct also for inverse hierarchy. 

We should stress that we have considered so far the CP conserving case HPS ansatz, 
with same-CP-sign neutrino mass eigenvalues, 

mi, ma, m3 > . (40) 

However, for all other CP combinations the denominators in Eqs. fl25ti27p are larger such 
that the deviations from HPS mixing pattern become smaller and correspondingly relax the 
bound in Eq. (fT5l) . In particular for the case of opposite CP signs there is no bound, as seen 
in right panel in Fig. [3l 

V. SUMMARY AND DISCUSSION 

We have studied the stability of the HPS mixing ansatz that could arise from a flavor sym- 
metry valid at some high energy scale, against supersymmetric radiative corrections (RGE 
running and threshold effects). We have adopted a model-independent minimal supergravity 
framework where supersymmetry breaking is universal and flavor-blind at unification. In 
this case we have found that the solar mixing angle can only be increased with respect to 
the HPS reference value. Under the assumption that all neutrino masses have the same 
CP-sign, this sets a rough upper bound on the mass of the lightest neutrino which, in turn, 
implies that the neutrinoless double beta decay rate is also bounded as a function of the 
mSUGRA parameters. In contrast, in the case of opposite CP signs there is no bound on 
the lightest neutrino mass. We have also shown that the atmospheric and reactor mixing 
angles remain essentially stable in all cases. It should not be surprising that the effect of 
radiative corrections is more important for the solar angle than for the others. It simply 
reflects the fact that the solar is the smallest of the two neutrino mass splittings. 
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We stress that in our approach we have assumed only that the matrix m^J'^'^ is diagonahzed 
by the HPS matrix at the unification scale and this gets modified only by minimal supergrav- 
ity radiative corrections, universal and fiavor-blind at unification. This concerns the struc- 



ture of the dimension- five operator, Fig. [H Additional radiative corrections 12 1 to the solar 
angle HPS prediction are expected, if the neutrino mass arises a la seesaw 
Their magnitude will be determined by the strength of the Yukawa coupling characterizing 
the Dirac neutrino mass entry in the seesaw mass matrix 19|. This will depend strongly 
on the details of the model, in particular, on whether Higgs triplets are present in the see- 
saw \V7\ or on whether the seesaw is extended [2^. Scrutinizing the schemes for which it 
is possible to decrease the solar mixing angle value predicted by the HPS mixing pattern 
towards its currently preferred best fit point value will be considered elsewhere 2l| , together 
with the related issue of the lepton fiavor violating processes that would be expected in these 
schemes. 
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